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active neutrino fluxes into the bulk, caused by the destructive interference 
from the Kaluza-Klein states. We discuss a class of models to illustrate this 
general feature. 

I. INTRODUCTION 

This article is dedicated to Prof. Hiroshi Ezawa on the occasion of his seventieth birthday. 
I met Hiroshi in the early 1960's, at the University of Maryland. Right away it was clear 
that I could learn much from him, both in physics and in mathematics. What I did not 
realize until later was his administrative talent and his superb quality in leadership, both 
amply demonstrated in his illustrious career. I would like to take this opportunity to wish 
Hiroshi a happy birthday, a happy retirement, and many many happy returns. 

The possible existence of extra (spatial) dimensions beyond our three was first suggested 
by Theodor Kaluza in 1919, later modified by Oskar Klein in 1926. When superstring came 
along, consistency requires it to live in six extra dimensions. Unfortunately, there is no 
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experimental evidence to date for the presence of these extra dimensions. That may be due 
to the smallness of the extra dimensions, too small even for the largest accelerators to see. 
The lack of Kaluza-Klein excited states up to about 1 TeV places an upper bound on the 
size of the extra dimensions to be about 10~^^ m. 

Inspired by the discovery of higher-dimensional D-branes in non-perturbative string the- 
ories [1], where open strings are trapped, a braneworld scenario was proposed [2] in which 
the Standard-Model (SM) particles are confined to our three dimensional world, called a 3- 
brane. Only SM singlets such as gravitons and right-handed neutrinos may leave our world 
to roam in the extra-dimensional bulk. In that scenario, SM particles have no Kaluza-Klein 
(KK) excited states simply because they cannot get into the bulk, so the upper bound of 
10~^^ m placed on the size of the extra dimensions is no longer valid. In order to probe 
the presence and the property of the extra dimensions, we must use either gravity, or the 
right-handed neutrinos i/r. 

A deviation from the inverse-square law of gravity can be used to detect the presence of 
an extra dimension. In a 3 -|- n dimensional world, the surface area of a sphere of radius r is 
proportional to r^"'"". Hence the gravitational force between two masses mi and m2 is given 
by G„m 1 777.2 /r^"*"", where G„ is the Newtonian gravitational constant in n extra dimensions. 
This is the case when r is less than the size R of the extra dimensions. Otherwise, the sphere 
is squashed in the extra dimensions to a size R, so the surface area of the squashed object is 
now proportional to r^i?". The resulting gravitational force G„(«;/i?")/r^ once again obeys 
the inverse-square law, where k is a computable geometrical factor. If a deviation from the 
inverse-square law is detected experimentally at r < R, then R marks the size of the extra 
dimension. Present experiments found no deviation down to about 0.2 mm [3], so R must 
be smaller than that. However, it could be as large as 0.1 mm. If so, and if there are at least 
two extra dimensions, then G„ = GgRn/ n is large enough for strong gravitational effects to 
be seen at TeV energies. This possibihty led to a lot of excitement and many papers. 

What if the size of the extra dimension is much smaller than 0.1 mm, or, there is only 
one extra dimension with such a large size? In that case gravity remains weak, and it is 
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powerless to yield any information in the near future on extra dimensions. 

What about the other probe, the right-handed neutrinos i^r? They are assumed to be 
absent in the SM. In any case, they are SM singlets, hence sterile, in the sense that they 
experience none of the SM forces. How can they be detected even when they are present? 

The answer is 'mass'. If vrS exist, they would probably produce neutrino masses through 
the Dirac-mass coupling iTrI'l with the left-handed neutrinos ul. A good indication of the 
presence of a right-handed neutrino is therefore the presence of a neutrino mass. 

No mass has been detected in the tritium /3-decay experiments. This places an upper 
bound of 2.2 eV on the mass [4] of the electron- ant ineutrino, Vg. However, the pioneering 
neutrino experiments led by Davis and by Koshiba, the beautiful data coming from Super- 
Kamiokande, SNO, and KamLAND [5], support an oscillation explanation for the missing 
neutrinos they detected. It demands at least two of three neutrinos to have non-zero masses. 
Specifically, if Mi, M2, M3 are the masses of the three mass eigenstates, then AM© = M| — 
Mf ~ (7.5 X 10-=^ eV)^ and AM^^ = |M| - M|| ~ (50 x 10"=^ eV)^. There are also 
astrophysical evidence suggesting the neutrino masses to be bounded above by 0.23 eV, if 
they are degenerate [6]. 

Unlike the quarks, oscillation experiments discover that neutrinos mix strongly among 
themselves. Three rotation angles and one phase angle are needed to describe the mixing of 
three left-handed fermions. In the case of quarks, all three rotation angles are small. In the 
case of neutrinos, two of them (^12 and ^23) are large and one of them (^13) is small. 

Now that wc know the neutrinos have a mass, we shall assume the right-handed neutrinos 
to exist. Where do we find them? Neutrino mixing is large but quark mixing is small; 
neutrino masses are small but quark masses are at least a million times larger. These 
differences suggest that neutrinos are quite different from the quarks. Since the left-handed 
quarks and the left-handed neutrinos behave in much the same way under the SM, that 
indicates the right-handed quarks and the right-handed neutrinos are very different. Right- 
handed quarks are found at the SM energies, this difference may be telling us that the 
right-handed neutrinos should be found elsewhere. 
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Where? The popular scenario is to assume the right-handed neutrinos to hve at very 
high energies. Through the seesaw mechanism [7], this scenario explains the smallness of the 
neutrino mass, though extra assumptions are needed to explain the large mixing of neutrinos 
this way. 

With the braneworld scenario, there is another possibility. Quarks are confined to the 
3-brane we live in, but urs are free to roam in the bulk. That distinction might give rise to 
the qualitative difference between quarks and neutrinos. In the rest of this article, we shall 
examine this possibility more closely. 

Each neutrino in the bulk yields an infinite tower of KK neutrinos in four dimensions, 
which we shall refer to as the bulk neutrinos, or the bulk states. These neutrinos are 
non-chiral, containing both the left-handed and the right-handed components. Since they 
originate from the bulk, they are sterile. 

To proceed further, let us keep two general questions in mind. First, can the braneworld 
scenario explain the difference between quarks and neutrinos in a natural way? That is 
important because that is the raison d'etre for going into extra dimensions. Second, do the 
data on neutrino oscillations even allow the extra-dimensions to exist? The second question 
is relevant because solar and atmospheric neutrino data demand the mixing with sterile 
neutrinos to be small, but in the braneworld scenario all the bulk neutrinos are necessarily 
sterile. 

The answer to these questions depends to some extent on the size R of the extra dimen- 
sion. In most of the recent literature [8], the size is assumed to be large. A size of i? = 0.1 
mm corresponds to a characteristic energy of 2 x 10"^ eV, putting it in the right ball park 
of the neutrino masses. Using perturbation theory, a weak coupling between the left-handed 
brane (the SM) neutrinos and the right-handed bulk neutrinos can be shown to produce a 
small neutrino mass as well as a small mixing with the sterile neutrinos. It however does 
not explain the large neutrino mixing in a natural way, though that can be arranged. 

What if the size of the extra dimensions is much smaller than 0.1 mm? Then the weak- 
coupling assumption explains nothing, so neutrino oscillation is no longer a useful probe 
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for extra dimensions. Neither is gravity in that case. We are then back to the unenviable 
position of having no way to tell the presence or absence of an extra dimension. 

That however is based on the assumption that the brane-bulk coupling is weak. We shall 
argue that it is more natural to expect that coupling to be strong, not weak. In that case, 
things are quite different, and neutrino oscillation can provide useful information even when 
the extra dimensions are small. 

Strong coupling with the bulk may induce a large mixing between the brane neutrinos, 
thereby explaining naturally why neutrino mixing can be large while quark mixing is small. 
This answers the first question posted before in the affirmative. But why then is one of 
the three neutrino mixing angles small, instead of being all large? We shall show that 
the smallness of that particular mixing angle is intimately related to the smallness of the 
mass-gap ratio AMq/AM^^^. 

What about the second question? With a strong couphng surely we expect a large 
amount of sterile neutrinos showing up in the the solar and atmospheric data, contrary to 
observation. Is there any way out of this fatal problem? 

There is. The details will be discussed in the rest of this article, but let us summarize 
here what is involved. 

The trick is to introduce an extra sterile neutrino on the brane, so that most of the 
mixings are between the sterile brane and bulk neutrinos, and not between the active and 
sterile neutrinos. In the strong coupling limit, the introduction of this sterile brane neutrino 
is forced on us by the dynamics; it is not arbitrary and artificial. This is so because one 
of the brane neutrinos is always absorbed by the KK tower of bulk neutrinos in the strong 
coupling limit. In other words, if we start with / flavor neutrinos in the brane, in the strong 
coupling limit there are only / — 1 mass eigenstates left on the brane. 

Therefore, in order to have two separate mass gaps needed to explain the solar and 
atmospheric neutrino experiments, we need to have three mass eigenstates in the brane, and 
hence four favor states to start out with. We know from the width that there are only 
three active neutrinos, Ve-i^in and so the fourth one must be sterile. We shall denote it 
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by i^s- 

With the extra dimensions small, the small mass of these neutrinos can no longer be 
explained in the usual way [8], so seesaw or some other mechanism must be invoked. In 
what follows we shall not ask the origin of these small masses, we simply introduce four 
parameters nia to describe the Majorana masses of the flavor neutrinos on the brane. 

All masses from now on are understood to be measured in some unit U, so the parameters 
nia and later on da are dimensionless. 

In addition to the four brane neutrinos, a minimal model would consist of a single 
massless bulk neutrino in five spacetime dimensions. It decomposes into an infinite tower 
of four-dimensional KK neutrinos, to be called the bulk neutrinos, or bulk states. With U 
properly chosen, the spectrum of this infinite tower can be taken to be the set of all integers. 

A Dirac mass coupling of the type dTT^VL + h.c. is introduced to couple the left-handed 
brane neutrinos ul to the right-handed bulk neutrinos ur. Since the KK tower comes from 
a single neutrino in the bulk, there is only one coupling constant da per brane neutrino. 

If the Dirac masses da are comparable to the Dirac mass of any of the charged fermions, 
and if the Majorana masses are comparable to the neutrino masses, then da is larger than 
mfe by more than a million times. Defining the overall coupling strength to be d^ — Y^^=i d'i-, 
and letting = da/d, it is therefore likely that Nature is operating in the strong-coupling 
regime, where d ^ m^, ea, 1. 

In summary, there are eight real parameters in the minimal theory. Four m^'s, and 
four daS. The fiavor neutrinos are the states when d = 0; the mass eigenstates in the 
strong-coupling limit are the neutrinos when d — > oo. 

To get the mass eigenstates for d 7^ 0, we have to diagonalize an infinite dimensional 
matrix, whose rows and columns are labeled by the four flavor neutrinos on the brane, and 
the inflnite number of flavor neutrinos of the bulk. The neat thing is that the eigenvalues 
and the eigenvectors of this inflnite dimensional matrix take on a very simple form in the 
strong coupling limit. 

As mentioned before, one brane neutrino is absorbed into the KK tower of bulk neutrinos. 
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Taking this into account, the final mass spectrum in the strong couphng hmit is as follows. 
The whole bulk spectrum is rigidly shifted by half a unit, so that the masses are now half 
integers. The three mass eigenvalues Mi,M2,M3 on the brane are sandwiched between 
the four Major ana masses rria of the flavor neutrinos. Namely, if we order the parameters 
according to mi < m2 < ms < 777,4, and Mi < M2 < M3, then ?77i < Mi < 7772 < M2 < 
rris < M3 < 772.4. We shall refer to this inequality as the ordering relation. It is a crucial 
feature of the strong coupling model. 

In the strong coupling limit, we are left with seven free parameters, four rUaS and three 
independent Cq's (because J2t,=i^l = It turns out that we can replace them by the four 
777a's and the three Mj's, provided the ordering relation is maintained. Technically this is 
quite important because it is much simpler to deal with the latter set of parameters than 
the former set. 

The eigenvalues, and hence the unitary overlapping matrix between the flavor and the 
mass eigenstates, can also be worked out. 

One can then compute the probability amplitude of a flavor neutrino oscillating into 
a flavor neutrino i/f,, after traversing a distance L. To do so, we must decompose the flavor 
neutrino Va into a linear combination of the eigenstates, because it is these normal states 
that propagate with a definite frequency. After a distance L, the mass eigenstates must all 
be converted back into the fiavor state to get the probability amplitude. 

When d is large, the fiavor states on the brane have only a tiny overlap with each of the 
bulk eigenstates. Nevertheless, since there are an infinite number of bulk eigenstates that 
a strongly coupled fiavor brane neutrino can reach into, the total effect of the bulk is not 
negligible. The contribution from the infinite number of bulk states destructively interfere 
with one another, so completely that any active neutrino that oscillates into the bulk will 
not be able to come back. In other words, the bulk acts like an absorber to the active 
neutrinos in the brane. 

Oscillations that go through the three brane eigenstates act just like ordinary oscillations 
without the presence of extra dimensions. 
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In other words, the effect of the extra dimensions is to cause part of the oscillating flux 
of active neutrinos to be lost in the bulk. This then is the signature of the presence of an 
extra dimension that we should look for. 

We can now describe the physical significance of the seven parameters in the theory. 
Mi,M2,M3 are the eigenmasses of the three active brane neutrinos. Since the mass is 
actually Mill, it is only these products that can be determined experimentally. Different 
values of U corresponds to different size R of the extra dimension, because U was chosen 
to make the flavor mass of the bulk neutrinos to be an integer. Since U alone cannot be 
determined from the experiment, neither can R. 

This is not to say that we can use this method to probe an extra dimension no matter 
how small it is. The strong-coupling requirement d ^ 1 places a limit how small R ~ 
can be, if we assume the Dirac mass dU ~ d/R to be comparable to the Dirac mass of 
the charged fermions. The precise value of course depends on what we use for d. Let us 
illustrate it with two extremes. If d/R is the electron mass 0.5 MeV, then we need to have 
3> 4 X 10~^^ m. If it is the top quark mass 175 GeV, then we can go down to an ^ 10~^^ 

m. 

The other four parameters, mi, 777.2,777.3,777.4 can be used to fit the three mixing angles, 
and the amount of absorption by the bulk. Note that there is only one free parameter to 
describe the potential absorption for any oscillating into any i/^, so there are predictions 
that can be potentially falsified. 

No absorption has been detected in the present data. Refined and precise data in the 
future may. It that happens, it is a good indication that an extra dimension exists. 

In the minimal medel, the absence of absorption can be achieved by letting 7774 — > 00. 
In that limit the mixing between the active neutrinos v^^n) and the sterile neutrino 
also disappears. Then we revert to a situation indistinguishable from the case without extra 
dimensions. In principle, there is actually a way that might tell them apart, because the 
three parameters 1711,1712,1713 in the minimal model are constrained by the neutrino masses 
Ml and M2 through the ordering relation. As such we may not be able to use them to fit 
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the three experimentally measured mixing angles. In reality these constraints are fulfilled 
in the fit, so these two cases become indistinguishable. 

When Ml approaches M2, the ordering constraint requires Mi = m2 = M2. This pinching 
of 1712 imphes ^13 = 0. So in this model, one obtains the interesting prediction that AMq = 
imphes ^13 = 0. The smallness of ^13 is then related to the smallness of the mass gap ratio 



This minimal model in five spacetime dimensions can be substantially generalized without 
changing any of the crucial features discussed above. 

These descriptions and conclusions will be put into mathematical formulas in the next 
few sections. In Sec. 2, the mass matrix, its eigenvalue, and its eigenvectors of the minimal 
model are examined. They are then used to compute the oscillating amplitude in Sec. 3. 
Generalization beyond the minimal model will be discussed in Sec. 4. 



The mathematical solution of the minimal model will be sketched here. For more details, 
please consult Refs. [9] and [10]. 

This model contains four fiavor neutrinos in the brane, and a single massless fiavor 
neutrino in a five dimensional spacetime. The latter decomposes into a KK tower of bulk 
neutrinos with integer masses, and the former are each given a Majorana mass TJia (a = 
1, 2, 3, 4). The brane neutrinos are coupled to the bulk neutrinos by a Dirac-mass coupling, 
with coupling strengths da- We assume all masses to be expressed in some common unit U, 
so that the eight real parameters rria and da arc dimcnsionlcss. Direct coupling between the 
brane neutrinos is assumed to be absent, and CP violation is ignored in this simple model. 

The symmetric mass matrix of this model is 



where m = diag(mi, m2, ^3, 7714) is the 4x4 mass matrix of the brane neutrinos, and 
B — diag(0, +1, —1, +2, —2, +3, —3, • • •) is the infinite dimensional mass matrix of the bulk 



AMI/AM^ 



II. THE MINIMAL MODEL AND ITS SOLUTION 




(1) 
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neutrinos. The coupling between the two is supphed by D, a 4 x oo matrix in which every 
element of the ith row is equal to 

The eigenvalue equation for the mass matrix is 

where w is a 4- dimensional column vector with components Wa, and v is an oo-dimensional 
column vector with components Vn- A is the mass eigenvalue. In component form, (2) reads 

rUaWa + daA = XWa, (3) 
b + {Bv)n = XVn, (4) 

where 

n 
4 

^ = H daWa- (5) 
a=l 

We shall choose the normalization of the eigenvectors by setting 6 = 1. 
The eigenvector components can be solved from (3) and (4) to be 
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A — n' 

Wa^A-^^{Ad)-^, (6) 
A - ma X-rUa 



where 



a=l 



i = E4 (7) 



a=l 

The constant A may now be computed to be 



^ " tr A - m tan(7rA) ^ ^ 
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The eigenvalue equation is obtained from (8) and (6) and the normahzation condition 6=1 
to be 

-,2 



a=l 



- tan(7rA) = V = dFriX) . (9) 

Let us solve this equation for the flavor eigenvalue (d = 0), and for the mass eigenvalue 
in the strong couphng limit {d — > oo). For d = 0, it follows from (9) that tan(7rA) = 0, which 
implies A e Z, unless A = rUa for some a. These are the expected eigenvalues because they 
are simply the matrix elements of the diagonal matrix M. when d—Q. 

For d — > oo, we should have tan(7rA) = oo, which implies A = Z + |, unless r(A) — 0. 
The former are the bulk eigenvalues, and the solutions of the latter are the brane eigenvalues 
Mj [i — 1, 2, 3). Since r(A) approaches ±oo when A — > ma + 0=*=, there is one zero of r(A) 
between each successive pairs of ma's. In other words, the ordering relation mi < Mi < 
m2 < M2 < m3 < M3 < m4 mentioned in the Introduction is obeyed. 

If d is large but not infinite, the eigenvalues will shift somewhat, but they are still 
bounded between consecutive ma's or consecutive n's. 

We will now show that the three independent parameters may be replaced by the 
three independent parameters , provided the ordering relation holds. The argument is 
based on the simple observation that r(A) is a meromorphic function of A, with four simple 
poles occuring at A = nia, and three zeros occuring at A = Mj. Moreover, r(A) approaches 
1/A when |A| 00. Hence we can write r(A) = IlLil'^ ~ ^a) / Y{t.=i{^ ~ ^o)- The residue 
at A = mft is then 

3 

e'^HK-MO/nK-^a)- (10) 

i=l a^b 

This formula determines once nia and Mj are known. To keep > 0, the ordering relation 
has to be obeyed. 

Let U\ be the normalized eigenvector, with components Ua\ = Wa/N and Unx = Vn/N. 
The norm N'^ of the original eigenvector {wa, v„) is given by A^^ = {Ad^s + T, where 
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^1 (A-m„)2' 



(A - nf 



(11) 



III. OSCILLATION AMPLITUDE 

Using (6) and (11), we can calculate the transition amplitude Aab from a brane neutrino 
of flavor h and energy E (measured in units of t/), to a brane neutrino of flavor a after it 
has traversed a distance L — 2Et (measured in units of U~^). The transition amphtude is 
determined by the formula 

Aabir) = Y: U:,U,,e-'''^ ^ A!,{r) + A^,{r), (12) 

X 

where is the contribution from the brane eigenvalues A = Mi, M2, M^, and A^ is the 
contribution from the bulk eigenvalues A G Z + |. 

When d ^ 00, the quantities Vn, Wa/Ad, s and T are all of order 1, so the magnitude of 
Wa is determined by Ad and the magnitude of is determined by {Ady. According to (9), 
Ad = l/{dr). For bulk eigenvalues, r = 0(1), so Ad = 0{l/d). This implies N"^ ~ T and 
Uax — 0{l/d). In that case the bulk components of an eigenvector are much larger than the 
brane components. For brane eigenvalues, A — 0(1), hence Ad — 0{d) and Wa — 0{d). In 
that case the brane components of an eigenvector dominate and iV^ ~ {AdYs. 

Let us denote the large- 0? value of UaUi by \4j, the value of s at A = Mj by s^, and 
1 1 {Mi - ma) by Xai- Then 

Vai = eaXail^i (1 < O < 4, 1 < i < 3), (13) 

and 

Al,{T) = jzv:yue-''^^'. (14) 

As it stands, ^ is a 4 x 3 matrix, but we can make it into a square 4x4 matrix by letting 
the last column to be \4/ = Co. The meaning of this last column will be discussed later. 
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Note that we can write T44 in the same form as the other Vai, namely, T44 = GaXaA/ -y/sl, 
provided we let A = 00. 

The resulting 4x4 matrix V can be shown to be real orthogonal. It depends on the 
seven parameters, nia and Mj. 

We have shown in (10) how to express in terms of these parameters. Similarly, it can 
be shown that 

= - 1[{M, - Mfe)/ n(M, - mj. (15) 

kj^i c=l 

We turn to the contribution from the bulk eigenvalues. Since Uax is unitary, it follows 
from (12) that Aab{^) = Sab, hence 

A^6(0)= 5ab-Ai{0). (16) 
Using (14) and the unitarity of the matrix V, we conclude that 

AabiO) = Sab - V:fVkf = Sab ' CaCf. (17) 

Therefore 

A^{0) = e„e,. (18) 

The contribution from the bulk eigenvalues can also be obtained directly from (12) and 
the paragraph following that equation to be 

^abir)= E jA^n '-Wx -,e-'''^^eae,F{T). (19) 

Ag^i {drrT (A - m„)(A - rrib) 

Both r and T are of order 1 as d ^ oo, so the contribution from each bulk eigenvalue to the 
sum is 0(1/0?^). Since there are an infinite number of bulk eigenvalues, the total contribution 
to the sum in (19) is not necessarily zero. In fact, we know from (18) that F(0) = 1 even at 
an infinite d. 

It can be shown that F(r) = g{K'^T), where K'^ = + tt'^cP), and that 
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is zero ai x — 0, and decrease to like (1 — i)/\^27rx for large x. This means that ^^(r) = 
eaCbgiK^r) is zero whenever r > 0, in the limit d — > oo. This function describes the 
absorption of the active neutrino flux into the bulk. 

Therefore, in the strong coupling limit, we end up with 

A.(r) = A!,{t) = E K*H.e-^^ \ (21) 
1=1 

In the limit — > oo, it follows from (10) that 64 — > 1, and hence from (7) that e;, — > 
for b = 1,2,3. Since the matrix V is real orthogonal, and 64 = V44, it also follows that 
Va4 = for a = 1, 2, 3. As a result, the active neutrinos I'c^n, do not mix with the sterile 
neutrino i/g, and the active neutrinos do not get absorbed by the bulk. 

If furthermore Mi — M2, then the ordering relation forces m2 — Mi. In that case 
V23 = 0. This means the mixing angle ^13 = if we identify the a = 2 flavor neutrino with 
Ue- Hence a vanishing AM|,/AM^^^ implies a vanishing ^13. 

IV. GENERALIZATION OF THE MINIMAL MODEL 

The final result (21) remains valid for almost all B in (1). This is reasonable because 
(21) does not depend on the property of the absorptive bulk, which B affects. For a detailed 
argument, please consult Ref. [11]. 

This research is supported by NSERC and FRNT. 
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